SEMI - RIEMANNIAN GEOMETRY WITH NONHOLONOMIC 

CONSTRAINTS 



ANNA KOROLKO, IRINA MARKINA 

Abstract. In the present article the geometry of semi-Riemannian manifolds with nonholo- 
nomic constraints is studied. These manifolds can be considered as analogues to the sub- 
Riemannian manifolds, where the positively definite metric is substituted by a nondegenerate 
metric. To study properties of the exponential map the Christoffel symbols and other differ- 
ential operators were introduced. We study solutions of the Hamiltonian system and their 
projections into the underlying manifold. The explicit formulae were found for a specific 
example of a semi-Riemannian manifold with nonholonomic constraints. 



1. Introduction 

Sub-Riemannian manifolds and the geometry introduced by bracket generating distributions 
of smoothly varying m-dimensional planes is widely studied interesting subject, which has 
applications in control theory, quantum physics, C-R geometry, the theory of principal bundles, 
and other areas. The main difference of the sub-Riemannian manifold from a Riemannian one 
is the presence of a smooth subbundle of the tangent bundle, generating the entire tangent 
bundle by means of the commutators of vector fields. The subbundle, which is often called 
horizontal, is equipped with a positively definite metric that leads to the triple: manifold, 
horizontal subbundle, and Riemannian metric on the horizontal subbundle, which is called 
a sub-Riemannian manifold. The foundation of the sub-Riemannian geometry can be found 
in \12\ [13| [HI 117] . The following question can be asked. What kind of geometrical features 
will have the mentioned triplet if we change the positively definite metric on the subbundle to 
an indefinite nondegenerate metric? We use the term semi-Riemannian to emphasis that the 
considered metric is nondegenerate in contrast to the positively definite metric, that referred 
as Riemannian. As it is known to the authors the present work is the first attempt to study 
systematically the geometry of semi-Riemannian manifolds with nonholonomic constraints, 
that we called the sub- semi-Riemannian manifolds or shortly ss-manifolds. In the present 
paper we study the exponential map and solutions to the Hamiltonian system that has no 
established terminology in the literature and can be called geodesies or extremals, see, for 
instance |12tll7j. The ss-manifolds have their own peculiarity that distinguishes them from the 
sub-Riemannian and semi-Riemannian manifolds. The simplest example of a semi-Riemannian 
manifold with nonholonomic constraints is the Heisenberg group equipped with the Lorentzian 
metric and it has been considered in [6l [TT] . It was shown in [11] that in contrast with the 
Heisenberg group with positively definite metric the Lorentzian type of the Heisenberg group 
possesses the uniqueness of extremals both of timelike or spacelike type. The structure of the 
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article is the following. Section [2] is devoted to main definitions. The collection of technical 
lemmas concerning Christoffel symbols is proved in Section [3j In Section [4] the extremals and 
exponential map are introduced, the extendability of extremals and Gauss lemma are shown. 
Some properties of the length are also studied. Section[5]is devoted to the differential properties 
of the exponential map. It is shown that the exponential map possesses an analogue of "local 
diffeomorphism" property, although is it not a diffeomorphism at the origin. The last Section [6] 
consists of the example of ss-manifold, where the explicit formulae of extremals are found. 

2. Main definitions 

Let M be a connected n-dimensional, n > 3, C°°-manifold. Let T x and T* denote the 
tangent and cotangent spaces at a point x G M, and (Y,£) the pairing between them, Y G T x , 
£ G T*. The tangent and cotangent bundles are denoted by T and T* respectively. Fix an 
integer m, 1 < m < n. Let S be a fixed subbundle of the tangent bundle T, S = \J S x , S x 

be a fiber over x, of the rank m. A subbundle S will be called bracket generating or complete 
nonholonomic, if the vector fields which are sections of S, together with all brackets span T x 
at each x G M. In this case any two points in M can be connected by a piecewise smooth 
curve 7(s) such that the tangent vector j(s) belongs to Sy/ 8 \ at each point j(s) where the 
tangent vector exists. The bracket generating subbundle S is called the horizontal bundle 
or horizontal distribution and a curve j(s) satisfying j(s) G SjU) 1S called the horizontal 
curve. A result of Chow [4], see also |16j . guarantees the connectivity of M by a horizontal 
curve. The necessary and sufficient condition on connectivity by curves tangent to a given 
distribution of a smooth manifold can be found in |18j . We notice that the connectivity of 
a manifold by horizontal curves tangent to a given distribution depends only on properties 
of the distribution and not on any metric defined on it or on the tangent bundle. If Y G S, 
let S + [Y, S] denote the subbundle of T spanned by S and all the vector fields [Y, X] , where 
X varies over sections S. A fiber at a point x G M is written as S x + [^(x), S x ] G T x with 
Y{x) G S x . Similarly we define bracket(fc, Y) inductively by bracket(2,y) = S + [Y, S] and 
bracket(fc, Y) = S+ [bracket(A; — 1, Y), S]. More generally we set bracket(2, S) = S+ [S, S] and 
bracket(/c, S) = S + [bracket(& — 1, S), S]. A restriction of a bundle to x G M is denoted by 
writing the subscript x, for example: bracket(/c, Y(x)) = S x + [bracket(/c — l,Y(x)),S x ] G T x . 
We will say that a vector field Y G S is a A:-step generator if bracket(A;, Y(x)) = T x for all x G M. 
Similarly, S will be said to be fc-step bracket generating distribution if bracket(fe, S x ) = T x for 
every x. From now on we work with a distribution S possessing the bracket generating property. 
By analogy with the sub-Riemannian metric we give the following definition. 

Definition 2.1. Let M be a smooth manifold, 5 be a bracket generating subbundle of the 
tangent bundle T. A sub-semi-Riemannian metric Q on S is a smoothly varying in x non- 
degenerate quadratic form Q x on S x . We abbreviate the long and tedious name of sub-semi- 
Riemannian metric by the term ss-metric. We call the pair (S, Q) the sub-semi-Riemannian 
(ss-Riemannian) structure on M. 

We remind that the index v of a metric is the maximal dimension of the space V x <Z S x , where 
the form Q x is negatively definite. If v = 1 then we call the ss-metric the sub-Lorentzian metric 
following the tradition in semi-Riemannian geometry. The ss-metric with the index v = is 
just a sub-Riemannian metric. Given Q x , we may define a linear mapping g x : T* — > T x as 
follows: for given £ G T x , the linear mapping W — ► {W, £), W G S x can be represented uniquely 
as W — ► Q X (W, X) for some X G S x , then X is chosen to be g x £. The map is called a 
cometric and completely defined by the two following relations: 
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(i) image of T* under g x is S x , 

(ii) g x and Q x are related by the identity 

(2.1) Qx(W,g x Q = (W,t) for all W E S x . 

Lemma 2.2. If Q x is symmetric, nondegenerate and has index v, then g x is symmetric, 
nondegenerate and has index v . 

Proof. We understand the action of the cometric g on T* x T* — > R (omitting x) as following: 
fiK^VO = (qC:^) f° r an y two covectors £ and ^ from T* . 

Thus by definition of the cometric g we have g{ip,^) = {gip,£) = Q(gip,gQ, which equals to 
Q(g£,gip) = (g£,,ip) = g(€,ift) by the symmetry of the ss-metric Q. 

Now, having the nondegeneracy of Q we prove the nondegeneracy of g, that is if ip) = 
for any ip E T* then £ = 0. First of all, we notice that the pairing (Y, £) is not degenerate 
by (|2.ip . Then, taking arbitrary ip E T* and setting y = g^>, we obtain 

5(^,0 = (y,0 = for any Fe5. 

This implies that £ = by the nondegeneracy of pairing and we conclude that <? is nondegen- 
erate cometric. 

Let 5^!" denote the kernel of g x , and C T* be the subbundle with fibers S^~. Then 
g x : T*/S x — > is bijection. The relation (12. ip shows that the index of Q x and g x coincides 
for any x £ M and that g x degenerates on S x . □ 

Conversely, given a symmetric nondegenerate linear operator g x : T* — > T x with image *S*x, 
there is a unique nondegenerate quadratic form Q x satisfying (|2.ip . We write gl k , j, k = 1, . . . n 
for the symmetric matrix defining the cometric g x to emphasis that it is a tensor of covariant 
type and operates with covectors. The matrix g x k is never invertible. 

A differential manifold M with a chosen subbundle S of the tangent bundle and with a given 
nondegenerate ss-metric Q on S will be called the sub- semi- Riemannian manifold or shortly ss- 
manifold. If the index v of Q is 1, then we call the triplet (M, S, Q) a sub-Lorentzian manifold 
and in the case of v = we get the sub-Riemannian manifold widely studied in (9j [T2l [Til [T7] 
and numerous references therein. 

We present a couple of examples of ss-manifolds. 

Example 1. 

Let us consider the following example of sub-Lorentzian manifold that we call the Heisenberg 
group with sub-Lorentzian metric. This example was considered first in [6j [7] and was also 
studied in We remind that the Heisenberg group H 1 is the space R 3 furnished with the 
non-commutative law of multiplication 

(x,y,z)(x',y',z') = (x + x', y + y', z + z' + -{yx' - xy')). 

This gives the R 3 the structure of a non-abelian Lie group. The two-dimensional horizontal 
bundle S is given as a span of left invariant vector fields 

d 1 d „ d 1 d 

X = ^ _ + o^« - ' Y= a o x ^~' 

ox 2 oz ay 2 oz 

that can be found as the left action of the Lie group. There is only one nonvanishing com- 
mutator [X,Y] = Z = We suppose that the Lorentzian metric Q is defined on S by 
setting 

Q(X,X) = -l, Q(yy) = l, Q(X,Y) = 0. 
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Thus the triple (R , S, Q) is called the Heisenberg group with the Lorentzian metric, and to 
differ it from the classical case H 1 we use the notation H^. We say the classical case bearing 
in mind the manifold (R 3 , S, d) with a positively definite metric d on S. 
The quadratic nondegenerate symmetric form Q on S is of the form 



Q = {Qab} 



-1 o 
1 



Take the basis of Lie algebra associated with the Heisenberg group, considered as the Lie group, 
(X, Y, Z) € T. The dual basis of T* consists of the forms dx, dy, ui = dz — \{xdy — ydx). We 
wish to find the cometric g = g J . Let gdx = a\X + C12Y . Making use of (|2.ip for W replaced 
by X and Y, we deduce that gdx = —X. In the same way we get g dy = Y. Thus, the equality 
g(C,0 = QidC^dO calculated for the basic forms implies the values g = gi k for j, k = 1,2. 
The rest of the terms vanish because of g{dx,uo) 
go; = 0. Finally we get 



(gdx,iv) = 0, g(dy,uj) = {gdy,uj) = 0, and 
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Example 2. Consider the example of ss-manifold related to the notion of Heisenberg- type 
groups based on quaternions [2 El [TO]. The manifold M is R 7 . The vector fields 
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form the basis of four-dimensional horizontal distribution S. These vector fields come from 
the infinitesimal action of the noncommutative group law multiplication 



L( x>z )(x',z') 



fx, z) o <x', z') 



[x + x', z + z' + — Im(x * x')) 



for (x, z) and (x', z') from R 4 x R 3 . Here Im(x * x') is the imaginary part of the product x * x' 
of the conjugate quaternion x to x by another quaternion x' . See the details in Section [6j The 
distribution S is bracket generating due to the commutation relations 

[Xi,X 2 ] = —Zi, [Xi,Xs] = Z3, [Xi,X^\ = Z 2 , 
[X2,X3\ = Z 2 , [X 2 ^Xi] = —Z3, [^35X4] = —Zi, 

where Zp = (3 = 1, 2, 3 form a basis of the complement to S in the tangent bundle. 
We define the ss-metric Q on S by the matrix 

/-l 0\ 
0-100 

QaP ~ 10 

\ 1/ 
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The ss- metric Q has index 2. The corresponding cometric g^ k is obtained like in the Example 1, 
has index 2, and assumes the following form 
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More details about the manifold of the Example 2 the reader can find in Section [UJ 

We refer to [3] for an example of sub-Lorentzian manifold based on the Lie group different 
from the nilpotent group. 



3. Christoffel symbols 



Recall that S% denotes the kernel of g(x) and S 1 - C T*, S 1 - = U x& mS^. The space S% is 
ie annihi 
geometry. 



the annihilator of S x in T* . From now on, we use the summation convention of the differential 



Lemma 3.1. Let v be a covector from S^. Then (v,Y) = for all Y G S if and only if 
g jk v k = 0, j = 1, ...,n. 

Proof. Let (v,Y) = for any Y 6 S. Pick up an arbitrary Y € S, then there exists a form u 
such that Y = goo. Moreover, we can assume that u> G T*/S ,_L . Indeed, if uj = uj\ + v\ , where 
uj\ € T*/S ± and v\ G S , then for any X G S we get 

Q(X,gu) =Q{X,gu l ) + Q{X,gv l ) = Q{X,gu l ) + (v x ,X) = Q{X,g< J j l ). 

Thus 

= (v,Y) = (v,guj) = Q(gcu,gv) for all Y = gu G S gv = 0. 

Here we used the symmetry and the nondegeneracy of Q. 

Conversely, having gv = we derive = (v,gu>) for any uj G T*/S ,_l . Thus, (v,Y) = for 
any Y = guj. □ 

Lemma 3.2. (a) If v is a section of the annihilator S- 1 -, then 



Jk dv k 



® dxP dxP k 

(6) If x(t) is a curve in M and v(t) is such that v(t) is a section of over x(t), then 



jk, v d 9 ] (x) 



dxv ±Pvk 



for all t (here the dot denotes the t- derivative). 
(c) If v and w are sections ofS ± , then 



Qgjk 
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d 

Proof. To prove (a) one applies — - to the identity g jk (x)vk(x) = which defines the null- 
bundle. 

To prove (6) we take the derivative — of the identity g jk (x(t))vk(t) = 0. 
Finally, to prove (c) first we apply (a) to obtain 

ik dv k dg 3k 

g J — — = — - — Vk- 

oxP oxP 

Then we multiply both sides by w € S and, making use of the symmetry of g, we get 

dg jk jk dvk kj dv k 

— — VkWj = —q J — — Wj = —a J Wj— — = I) 

dxP 3 y dxP 3 y 3 dxP 

because of g kj Wj = 0. □ 

The following question can arise: how the information about the bracket generating prop- 
erties reflects in the properties of gl If X, Y € S, then [X, Y] is an element of T/S. If X £ S x 
then there exists £ € T* with X = g x £, and similarly Y = g x rj. The covectors £ and r\ are not 
defined uniquely, as it was shown in the proof of Lemma 13.11 Thus they should be regarded 
as elements of T*/S^. The annihilator S 1 - contains all the necessary information concerning 
the commutators through the pairing ([X, Y],v), when v varies over S- 1 -. Let us consider the 
trihnear form ([g^gr,],v) on (T*/S^) x (T*/S^) x Sjr. 

Lemma 3.3. In local coordinates 

( ■ da rq ■ da rp \ 
(3-1) U,gri],v) = [g»Jjg. - g>*JL.) t pVq v r 

for v varying over S x for any x. 

Proof. Let £ and rj denote any sections of T*. Then X r = <? rp £ p and Y r = g rq r\ q are sections 
of S, and 

[X ,Y Y = Xi JL r -Y,JL X * 

Taking the inner product with v € S^, we find that the last two terms are annihilated since 
g rq w^v r = g qr v r ^H = and g rp ?^v r = g pr v r ^k = 0. Thus we obtain □ 

OXl OX3 OXl Ox3 

We want do define the analogue of the Christoffel symbols but with the raised indexes and 
see the relation between them and the trilinear form defined in Lemma 13.31 We write 

For sections £ € T* and v 6 S 1 - define r(£, v) £ T by T k {^,v) = T k P q ^ p v q . In classical case 
of differential geometry the Christoffel symbols are used to express the covariant derivative in 
local coordinates. Unlike to the classical covariant derivative, which associates for two vector 
fields another vector field, the operator T, as we will see from the following lemma, associates 
a vector field for a pair of covector fields and, moreover, the resulting vector field is horizontal. 

Lemma 3.4. r(£,u) is a well-defined vector field; that is it is independent of the choice of 
coordinates. Moreover, T(£,v) is a horizontal vector field andT(£ + w,v) = T(£,v) forw £ S 1 - , 
so that V : (T* /S ± ) x S 1 -> S. 
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Proof. Let us prove that T k (^,v) transforms as a tangent vector at each x: T k (y) = T d (x)^^, 
where y = ip(x) and ip is a local diffeomorphism determining a new coordinate system. By 
g, £ and v denote the expressions for g, £ and v in the new coordinates. We have 

^ dx* dx k 9 [V) ~ 9 [X) dxP dx q ■ 

In the new coordinates 

~ w dg kq = ( ab dy^df_\d_( cd dy k dy q \ dx l 



dyi \ dx a dx b J dx l \ dx c dx d J dyi ' 
hence 

■ dg kq ~ „ [ ab <9y J 3y p <9g cd <9y fc ch/ 9 <9x z 



<9x a dx b dx l dx c dx d dyi 



(3.4) 



+gab dyidy^d^ g cd ( d 2 y k dy q + dy k 3 2 y q \ 
dx a dx b dyi \dx l dx c dx d dx c dx l dx d J 



The first term of (|3.4p equals to 

Qy k / . rliii f)n c d Rt^ \ / Rn^\ / , f)n°d\ 



dx c 



ab ^d^dx^\ = ( d]f \ f ab dg^\ 
dx a dx l dyi J \ dx c J \ dx a J 



Changing indexes b to p, d to q, and a to j we recognize the tangent bundle transformation of 
. dg kq 

g 3V -z — -£ v v a . The middle term vanishes since 

0X3 

(3.5) f d ^v q = f d v d = 

and the last term gives 

The middle term in T kpq £ p v g transforms as follows 



dyi ipVq -\dxc) \ 9 dx- J UVb -\dxc) \ 9 dxi I "■ 



The other terms vanish by the same reason as in (|3.5I) . The third term in T pq £ p v q in the new 
coordinates takes the form 



dy k ( ic dgP q \ lb dy k cdt d 2 y q 



(3 - 7) - w w^ q = ~d* ) ~ 9 dJ f ic dxW~ V - 

We see that the last term from (|3.7p is canceled with (|3.6j) (after the change of indexes). Taking 
together the rest of terms, we get the desired transformation law 



To show that r(£, v) 6 S we take a covector uj G S% and calculate 

{r&v),u) =T k P% )Vq uj k . 
Using ()3.2p and Lemma 13.21 we argue for each term of T k (^ p v q ) as it follows 
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and get (T(^,v),u) = 0, that implies r(£,u) G 5. 

The property T(£ + to, v) = v) for lo G S x follows from 

(3.8) T ^ v) = ^ ir ^ + ^£ ipVt 

and Lemma 13.21 □ 
Analogously to sub-Riemannian situation [T7] we have 



Theorem 3.1. A vector field X G S is a 2-step bracket generator if and only «/T(£, •) : S 1 - — > S 
is injective, where X = g£. In particular, S satisfies the 2-step bracket generating hypothesis 
if and only if •) : S 1 - — > S is injective for every nonzero form £ G T*/S ,_L . 

Proof. In the proof we exploit the properties of different linear mappings which we defined up 
to now. We have 

(3-9) (\9Z,gri\,v)=zr(£,v)v 

by (|3.1h and (|3.8j) . In order to show that the vector field X is a 2 step bracket generator we 
must show that the vector fields [X, Y] mod S fill out all T/S (at each x) as Y varies over S. 
In other words, the mapping 

(3.10) [X, ■] mod S : S — > T/S is surjective at each x. 

Since at any x the space T* / S 1 - is canonically isomorphic to the dual to S, statement (|3,10D 
is equivalent to 

(3.11) [<?£,<?(•)]: T*/S ± — > T/S 1 is surjective at each x, where X = g£. 

We notice that at any x the space S 1 - is canonically isomorphic to the dual to T/S. Thus (|3.1ip 
is equivalent to 

(3.12) T(£, •): S 1-1 — > 5 is injective at each x. 

□ 

We discussed earlier the relation between classical notion of covariant derivative and the 
Christoffel symbols. The closest notion to the notion of covariant derivative is symmetrized 
covariant derivative that was defined in |17j . It is natural to define the same concept on 
ss-manifolds. 

Definition 3.5. The symmetrized covariant derivative Vsym of a vector field Y is defined by 

, udY q dY k dq kq 

(3.13) (VvmY f, = g >,— +g „—- Y ,^. 

Thus (Vsym)*.: T X ^T X X 

Lemma 3.6. Vsym is a well-defined differential operator from tensors of type (1,0) to sym- 
metric tensors of type (2,0). Furthermore, ifY is a vector field from S, that is Y = g£, then 
(VsymY) kq v q = 2T k ^,v)foranyveS ± . 

Proof. The symmetry follows from the symmetry of the cometric g. Let us show that Vsym^i 
Y G T, transforms as a tensor field of rank (2,0). We check how the first term of (|3. 13|) 
transforms with the following change of coordinates: 

dxi ' dxP dxi ' 
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dyj dx a dx 13 dx l \ dx a J dyj 

a0 dy k dy j dx l d 2 y q ya | a0 dy k dy j dx 1 dY a dy q 
9 dx a dxl 3 dyi dx l dx a 9 dx a dx 13 dyj dx 1 dx a 
aidy k ay ya | al Qy k dY a dyQ 
dx a dx l dx a dx a dx 1 dx a 

Analogously, the second term 

wi^l = ^ dyq d2yk v a + </* dyQ dya dyk 



dyj dx a dx l dx a dx a dx 1 dx a 



And the third term 



yj d 9 kq _ dy J ya d f ab dy k dy q \ dx 1 



dyj dx a dx 1 \ dx a dx b J dyj 

_ dyj y a dg ab dy k dy q dx 1 
dx 01 dx 1 dx a dx b dyj 

, ^_ Y a ab dy9 9x1 | dyj Y a ab dyk d2y<} d °° l 

dx a dx l dx a dx b dyj dx a dx a dx l dx b dyj 

After summation and necessary renaming of indexes we have 



Vsym ' ~ da*dxj"' 8ym } ' 



\kq 



Let us show the second statement of the theorem. We assume Y = g^ € S and v 6 S . 
Then from (JOj) follows that 

(v sym ^) fc s = 0^&+^v*g>>, 

+ \ 9 dx^ p + 9 9 dxj) Vq 

Qgkq 

~dxT Vq 

= 9*%9 k ^ + 9 kj ^ P v q = 2r k (t,v) 
by d33J and gv = 0. □ 



9 j % 



4. HAMILTONIAN SYSTEM, EXPONENTIAL MAP AND LENGTHS OF CURVES 

The distribution S x at each point x of ss-manifold M has the structure of M. m equipped with 
a nondegenerate metric Q x of index v. The presence of the nondegenerate metric yields the 
following trichotomy 

Definition 4.1. A horizontal tangent vector w € S x is 

spacelike if Q x (w,w)>0 or w = 0, 

null if Q x (w,w)=0 and w ^ 0, 

timelike if Q x (w,w) < 0, 

nonspacelike if it is either timelike or null. 
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The set of all null vectors in S x is called null-cone at x G M. The category into which a 
given tangent vector falls is called its causal character. The terminology is adapted from the 
relativity theory, and particularly in the Lorentz case, null-vectors are called lightlike. For the 
nice and complete presentation of the semi-Riemannian geometry see |15j . 

The covectors G T*/S x receive the same causal structure according to the values of 
(g x ^(x), The covectors v G S' L (x) we shall call annihilators to distinguish them from the 

null-covectors. 

Definition 4.2. A horizontal tangent vector field X G S is spacelike, null or timelike if at 
each point x € M the vector X(x) is spacelike, null or timelike respectively. 

Definition 4.3. A section £ G T*/5 _L is spacelike, null or timelike if at each point x G M the 
covector is spacelike, null or timelike respectively. 

As we mentioned from the beginning, we work with the special class of admissible curves that 
tangent to the distribution S and that we called horizontal curves. We borrow this name from 
the sub-Riemannian geometry. We say that a horizontal curve c(s) is spacelike, null or timelike 
if the tangent vector c(s) is spacelike, null or timelike respectively at each point of c(s) where 
it exists. We can give the definition of the spacelike, null or timelike curve using the causal 
structure of the cotangent space T* according to the sign of (y c ( s )£(s), £(s)) = Q c ^(c(s), c(s)), 
where c(s) = <? c (s)£( s )- We call a horizontal curve the causal if the tangent vector c(s) (the 
covector £(s)) is nonspacelike 

In the sub-Lorentzian case we also introduce (as in the classical Lorentz manifolds) the time 
orientation. 

Definition 4.4. A time orientation on (M, S, Q) is a continuous horizontal timelike section T 
os S. 

If M admits a time orientation T, then T divides all nonspacelike horizontal vectors into 
two disjoint classes, called future directed and past directed. Namely, nonspacelike w G S x is 
said to be future (respectively past) directed if Q x (T(x),w) < (respectively Q x (T(x),w) > 0). 
We assume that any considered in the article sub-Lorentzian manifold (M, S, Q) will be time 
oriented. 

Since g: T* / S 1 - — > S is injective the time orientation can be brought to T*/S'- L . 

Definition 4.5. The globally defined section r G T* such that T = gr is time orientation on 
T*/S L . 

The covectors from S 1 - we can consider as null-covectors. 

The notion of arc length of a curve segment in Euclidean space generalizes in a natural way 
to ss-manifolds. Since the term "arc length" can be misleading since, for example, a null-curve 
has length zero. Therefore, we use the name "natural parameter" in stead of "arc length" . 

Definition 4.6. Let c : [a, b] — > M be a piecewise smooth curve segment in a ss-manifold 
(M, S, Q). The natural parameter of c(s) is 

L(c)= [ b \Q(c( S ),c(s))\ l / 2 ds. 

J a 

As in the classical case it can be shown that 

• the natural parameter is not changing under the monotone reparameterization and 

• if c(s) is a curve segment with |c(s)| = \Q(c(s), c(s))\ 1 ^ 2 > 0, there is a strictly increasing 
reparameterization function h such that 7 = c(h) has I7I = 1. 
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In the latter cases 7 is said to have unit speed or natural reparameterization. 

Now we define the extremal using the Hamilton function. Given the cometric g x : T* — > S x 
we form the Hamiltonian function 

(4-1) H(x,0 = ±{g x (0,0 

on T*. To emphasize the dependence of the cometric on x we write g{x) instead of g x when 
it is necessary. If we have the orthonormal basic X\, . . . , X u , . . . , X m of S we can write the 
Hamiltonian function in the form 

(4-2) jff(x , e) = __^(x j (x),0 2 + 2 S ( X j( x )>&> 

j=l j=v+l 

where v is the index of g x . Consider the Hamiltonian equations 

x(s) = V^H{x, 0, i(s) = -V x H(x, 
that explicitly can be expressed as 

x k ( s ) = g k ^x(s))^(s), k = l,...,n, 

(4-3) Us) = -l ^^h pisMs). 

An absolutely continuous curve T(s) on M satisfying (14.31) is called a characteristic of H . In this 
paper we will consider only the bicharacteristics T(s) such that H(F(s)) = H(x(s),£(s)) 7^ 
that are called in literature the normal biextremals. The detailed discussion of the structures 
of normal and abnormal geodesies see, for instance, |12[ W3[ [T|. Since we work only with 
normal biextremals we will drop the word "normal" for shortness. If H E C 1 (T*) then an 
extremal, of H is a curve x(s) which is a projection on manifold of some biextremal T(s) 
of H. The bicharacteristics of a Hamiltonian H € C k (T*) are curves of class C k along which 
H is constant. In this case it means that an extremal has a parametrization by the natural 
parameter. The next result is the consequence of this. 

Proposition 4.7. Ifj: [a,b] — > M is a normal extremal, then either Q^^( r y(s),j(s)) < or 
Q'r(s)( A f( s )^i( s )) = or Q 7 (s)(7(s),7(s)) > for all s € [a, b]. Moreover, if 7 is nonspacelike 
in the sub-Lorentzian manifold, then it does not change its orientation. 



Proof. We have 



2"w wv " Iy " 2 

which is constant along 7. The orientation preserving property of a smooth curve is obvious. 

□ 

It is possible to reformulate Proposition 14.71 in terms of cometric g. 

Proposition 4.8. Ifj. [a,b] — > M is a normal extremal, then either (£(s), <7-y( s )£( s )) < or 
(C(s),9 7 ( s )C(s)) = or (C(s),5 7 (s)C(s)} > for all s € [a, b}. 

Let us define the energy for the curve c : [a, b] — > M by 
(4.4) E[c) = ! \Q(c(s),c(s))\ds. 



In semi-Riemannian geometry extremals 7(5) are defined as curves which have parallel tangent 
vector field (7) or, equivalently, which have the acceleration zero: 7(5) = 0. It is true that semi- 
Riemannian extremals lift to solutions of (|4.3|) on the cotangent bundle. Thus the definition 
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of extremals like Hamilton extremals is correct generalization. For the sub-Riemannian and 
sub-Lorentzian cases see [811144 IT?]. Also, if we formulate the variational problem of minimizing 
energy E(c) over all smooth horizontal curves joining points p and q in M then the associated 
Euler equation is (|4.3p . Notice also that if we differentiate the first equation and substitute 
the second we obtain 

(4.5) ^ {s) + T k (t,0=0 

which is the analogue of the equation of the extremals in semi-Riemannian geometry. Notice, 
that we can not solve (14, 5h for £ in terms of x in any way. Thus (14. 5j) does not reduce to the 
equation in x alone. Neither (|4.5h together with x = g£ is equivalent to (|4.3j) . 

Given p € M, u € T*, and the coordinate system with the origin at p, the existence and 
uniqueness theorem for ordinary differential equations guarantees that the solution exists and 
is unique on an interval around zero provided the initial conditions x(0) = p, £(0) = u. As 
on sub-Riemannian manifolds the solution to (|4.3p can be continued as long as x(t) remains 
in M. 

Lemma 4.9. Let x(s) be a normal extremal for < s < a and suppose x(s) remains inside a 
compact subset of M. Then x(s) can be extended beyond s = a. 

Proof. Over the compact set K C M choose an orthonormal basis v^(x), . . . , v( n ~ k \x) of St- 
andi complete to an orthonormal basis of T* by adding u^(x), . . . ,u^ k \x). By definition all 
sections are smoothly varying on the compact set and hence bounded. Then the section £(s) 
along the extremal x(s) can be written as 

m n—m 

(4.6) £(*) =^ aj (s)u®(x(s)) + h(s)vW(x(s)), 

3=1 1=1 

where m is the rank of S. Consider (|4.3p as a system of equations for x(s), a,j(s), and bi(s). 
The functions x k (s) are uniformly bounded on K. Let us show that the functions Oj(s) and 
bi(s) are also bounded. We have 

1 m 

H(x(s),Z(s)) = -Y,(9^ ] ^)a 3 {s)a k {s) 

3=1 

by (|4.6p . Since extremals do not change the causal character and the Hamilton is constant 
along them, the value of the matrix g(u^'(x(s)),u^ k '(x(s))) is bounded from zero on K. It 
follows that aj(s) are uniformly bounded along extremals. Let us show that bi(s) are bounded. 
We write as 

m n—m 

(4.7) 600 =^2a J (s)u { i\x(s)) + £ 6,( a )««(x(.)), 

3=1 1=1 

where v/P and are coordinates of and in the local chart coordinates. We substi- 
tute (|4.7p in the second equation of (|4.3p and take into account the first one also. Notice, that 
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the terms involving gv and _ , v v v g vanish by Lemma [3 .21 c) since v G S ± . Finally, we get 

m n—m m o (j) n—m m a (J) 

E + E U w + E E + E E ">%9"a t uf 

j=l j=l j=l 1=1 j=l 1=1 

-. rj pn wn m m n—m n—m m 

-2%k(L E + E ^ E H w + E E 

i=i j=i i=i i=i j=i z=i 

If we dot both sides of equation with , uj^ , . . . , m ' , then we obtain the linear system 

b = Ab + C, 6 = (6i,...,& n -m), 

where the matrix A and the vector function C linearly depend on bounded functions a,-, 
and hence A and C are bounded. The linear system of the first order differential 
equations with bounded coefficients has bounded solution. We conclude that bj are bounded 
for j = 1, . . . ,n — m. Thus all the functions x k (t) and are uniformly bounded, and the 
local existence theorem implies the solution of (14. 3p extends. □ 

Now we can define the exponential map. 

Definition 4.10. If p G M, let D p be the set of covectors w in T* such that the extremal 
x w (s) is defined at least on [0, 1] and x(0) = p, £(0) = w. The exponential map of M at p is 
the function 

exp p : Dp — > M, such that exp p (w) = 

The set L> p is the largest subset of T* on which exp p can be defined. Fix w £ T* and r£l. 
Then the extremal s — > x w (st) is such that t£(0) = rio. Hence x ru ,(r) = x w (st) for all r and 
s where the both sides are well defined. Particularly 

exp p (sw) = x sw (l) = x w (s). 

As in the sub-Riemannian geometry the exponential map is always differentiable, since the 
solution of the Hamiltonian system depends smoothly on the initial data. But the exponential 
map is not a diffeomorphism at the origin. The reason is that all the extremals emanating 
from p must have tangent vectors in S p , but for any annihilator v G S 1 - we have 

exp p (v) = x v (l) = p, since x k = by (|4.3p . 

We prove the following analogue of Gauss lemma. In lemma we use the identification of a 
cotangent space T* at p with the tangent to T* space T U (T*) at point u G T*. The covector 
w G T p * at point u G T* is identified with the vector w G T U (T*). The radial vector r G T U (T*) 
means that it is a scalar multiple of a covector u G T*. 

Lemma 4.11. Let u be a cotangent vector in T* such that and lies inside D p . Let r be 

a radial vector and w be any other covector at point u G T* Then 

(i) 

(g p r,w) = Q exPpU (d(exp p ) u w : d(exp p ) u r) 
provided d(exp p ) u w G 5 ex p p u 

(ii) 

(g p r, w) = (d(exp p ) u w, £} 
where £ is a cotangent lift of the extremal t i— > exp p (tu) at t = 1. 
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Proof. Let us prove (i). Since r is radial, we can assume r = u. Take the curve v(s) = u + sw 
in T*. Let us suppose that the exponential mapping is defined in the cylindrical neighborhood 
Dp x [0, 1]. Consider the parameterized surface x : A — > M, A = {(i, s) : ^ t ^ 1, — e < s < e} 
given by x{t, s) := exp p (t(u + su>)). Note that 

d(t(u + sw)) d(t(u + sw)) 

Ws (1 ' 0) = w ' at (1 ' 0)=tt ' 

and the curves t i— > x(t,s) are extremals for any fixed s starting from the point x(0, s) = p 
with the initial covectors u + sw. Then 

dx dt(u + sw) 

— (1,0) = d(exp p ) u — (1,0) = d(exp p ) u w, 

dx dt(u + sw) 

— (1,0) = d(exp p ) u — (1,0) = d(exp p ) u u. 

dx dx 

Thus, we need to show (g p u,w) = Q C xn u(t^> -tt) (1, 0). Let s) be a cotangent lift of 

p os Ot 

dx dx 

the extremal t i— ► exp p (to(s)), particularly £(l,s) = it + sw. Then Q e xp p u(^ _ ) ■^■) (1, 0) = 
9x 

(—,0(1,0) and our main aim becomes to show 
os 

(4.8) (^,w) = (— ,0(1,0). 

We denote f{t,s) = {—,£)(£, s) and calculate the derivative — (t,0). We have 



3/„ \ 9& f . ,dx k d 2 xP 



Replacing ~~Qj~(t' s ) anc ^ ~~Qj~(t' s ) from the Hamilton- Jacobi equations (|4.3p we obtain 

a/^ = _ i gg!!^ g) gp(t> 8)gg(t> a) fo^ (ti a) + gp(t> ^^(gw^ft, 8) )^( f> 8)) 

(4.9) = j^(s(aOf (*,«),£(*,«)}) for any f and a. 

Since the Hamilton ^(g(x)^(t, s), £(t, s)) is constant along the extremal, then (|4.9p can be 
written as 

Then 

9/ 9/1 \ 

-^(^,0) = — (^-(5p(« + sw), (n + sw))J(t,0) = (c/ p n, w) for any t. 

We have 

9x 

/(0,0) = (—,0(0,0) = (d(exp p ) u tw,n}(0,0) =0 

that implies 

dx 

f(t,0)=t(g p u,w) /(l,0) = (— ,0(1,0) = {g p u,w). 

We proved f|4.8j) and hence (i). 
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To prove (ii) we argue in a similar way. Take the curve v(s) = u+sw in T* and parameterized 
surface x(t, s) = exp p (i(u + sw)). Let £(1, s) be a cotangent lift of the extremal t t— » exp p (tv(s)) 

dx 

at t = 1. Since d(exp ) u k; = —(1,0) the statement (ii) is reduced to (14.8ft . □ 
F OS 

Let c(t) be a C 1 piecewise curve in M for i € (a,b), where (a,b) is an interval in R. We 
remind that a curve c(i) is called horizontal if c(t) 6 for any t £ (a,b). A section £(t) is 
called a cotangent lift of c(t) if £(i) € an d fl^ = x{t) for every i where it is defined. The 
notion of the natural parameter or arc length (|4.6p for c(t) : (a, 6) — > M can be reformulated 
as follows 

b 

l( C ) = j(g c{t) m,m) 1/2 dt. 

a 

Let us focus for the moment on the case of sub-Lorentzian manifold. At each point p G M 
the distribution S p and the cotangent subbundle T* / Sp carry the structure of the Lorentz vec- 
tor space and thus inherit the typical features of the Lorentz structure. Since the orthogonal 
complement w 1 - to any timelike vector w is spacelike then the vector space S p can be decom- 
posed into the direct sum Km © w. The same regards the cotangent vector space T* /Sp. We 
define the future timecone in S p by 

C(T(p)) = {X(p) G S p : Q p (T(p), X(p)) < 0}, where T is the time orientation on S p . 

Analogously the future timecone in T* is 

C( r (p)) = { w £ T* : {g p w,r(p)) < 0}, where r is the time orientation on T* / S p . 

There is a consequence that vectors (covectors) v,w are timelike if and only if Q(v,w) < 
{{gv,w) < 0). In vector spaces with positively definite metric the Schwarz inequality permits 
the definition of the angle 9 between v and w as the unique number < 9 < tt. The analogues 
Lorentz result is as follows. 

Proposition 4.12. Let v and w be timelike vectors in a Lorentz vector space equipped with 
the scalar product (v , w) . Then 

(1) > where \v\ = | (v, t>) I 1 / 2 , \w\ = | (w, u;) I 1 / 2 . The equality is possible if 
and only if v and uu are collinear. 

(2) There is a unique number i? > 0, called hyperbolic angle between v and w, such that 

{v,w} = — \v\\w\ cosh 

Consider a piecewise smooth timelike curve c(t). The timelike means not only that every 
c(t) is timelike, but that at each break t{ of c 

Qc(u)(c(t;),c(4) <o. 

Here the first vector derives from c on the interval before break, and the second from 

the interval after break Thus c does not switch timecone at a break. Similarly we 

require that a piecewise smooth causal curve does not switch causal cones at a break. 

Lemma 4.13. Let p be a point at Lorentz manifold M. Suppose that 7: [0,b] £ T* is a 
piecewise smooth curve starting at the origin such that a = exp p o 7 is timelike. Then 7 
remains in a single timecone of T* . 
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Proof. We consider two cases. The first one is related with the smooth curve 7 and the second 
case will be general. Thus, we assume that j(t), and hence a(t) are smooth in all the domain 
of definition. Since ^(7(0), 7(0)) = (<? P 7(0), 7(0)} = Q p (a(0), d(0)) < 0, then the curve ^(t) is 
in the same timecone for t € (0, e), where e is sufficiently small. We also conclude that -y(i) 
maintains in the same timecone with t G (0, e) for sufficiently small e > 0. Let us denote by 
r 7 ( t ) the radial tangent vector in T^m(T*) corresponding to the timelike covector j(t). The 
vector r 7 ( t ) is timelike and therefore, (g p j(t), ?" 7 (t)) is negative for t G (0, s). We calculate 

(4-10) j t (g pl (t), 1 (t)) = 2{g p j(t), 1 (t)) 

is negative for t € (0,e). Since a(t) = d(exp p ) 1 ( t \ A f and it is in S exPp jM the last expression 
in (|4.10p is equal to 2Q exp 7 ( t ) [c\(t), d(exp p )y^) r j(t)) by Lemma [4.111 We conclude that so 
long as 7 remains in timecone the radial vector 77 and the vector d{exp p ) 1 r t \r 1 n\ remains 
timelike. Thus Q exppl ( t ){a(t),d(exp p ) l{t) r^ t) ) hence (g p j(t),~f(t)} and hence f t {g p l{t)^{t)) 
remains negative. But 7 can leave the timecone only by reaching null-cone or the origin. In 
any of these cases {g p "j(t), 7(t)) = 0. Thus 7 must remain in the same timecone. 

Now suppose that 7 and hence a is piecewise smooth. We know from the first part of the 
proof that on its first smooth segment 7 stays in the same timecone and therefore in the first 
break 

(g a (t )i(ti),r(5( tl )) < 0. 

Hence by Lemma 14.111 

Qeo;p P 7(t 1 )(«(*r)^( ea; Pp)7(ti) r 7(ii)) < °- 
The additional condition on a at breaks keeps a(tf ) in the same timecone, namely at t\ = 
d(exp p )r f u 1 )rr f Ui)- So, again by Lemma [4.111 (g a (t 1 )'y(ti), r p(ti)) < an d therefore it follows as 
above that 3^(57,7) can not change signs at breaks. Hence the argument for the smooth case 
remains valid. □ 

Minor changes in this proof show that the lemma remains true if the words timelike and 
timecone are replaced by causal and causal cone. 

Lemma 4.14. Let U be a normal neighborhood of p in a Lorentz manifold. If there exists a 
piecewise timelike curve a inU from p to q, then the segment a of the extremal from p to q is 
the unique longest timelike curve in U from p to q. 

Proof. We understand the uniqueness as the uniqueness up to a monotone reparametrization 
and can suppose that a is parameterized by arc length. If a : [0, b] — > U is a timelike curve 
in U from p = a(0) to q = a(b), then from the lemma above the lifting w(t): [0,6] — * T*, 
a(t) = exp p ow(t) remains in a single timecone in T*. The section w(t) is the timelike vector 
for any t £ [0,1] and therefore define a unite timelike section u(t) = \(g p w(t), w(t))|~ 1//2 u;(t). 
Since (g p u(t) , u(t)) = the vector field u(t) is spacelike. Let us write r(t) = |(g p zi;(t),u)(t))| 1 / 2 , 
then w{t) = r(t)u(t). We calculate 

a(t) = ^exp p w(t) = d{exp p ) w{ t)r(t)u{t) + d(exp p ) w ^r(t)u(t). 

Let us denote by £(t) the cotangent lift of the extremal 7 : s 1— > exp p (su(t)) at s = r(t). Since 
u(t) is timelike the same does £(t). Since a is horizontal, then d(exp p ) w u\r(t)u(t) G S a ( t y 
It means that there is spacelike section r)(t) such that d(exp p ) w ^r(t)u(t) = gai^t) 7 !^) almost 
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everywhere and orthogonal to the section £ by Lemma [4.11i Then the length of timelike vector 
a is given by 

\ 1/2 / \ 1/2 

Qa{t)(a,a)J = (f(t)£{t) +r],r(t)g a{t) £(t) + g a(t) r](t))j 

j(t)\ 2 ~ (9*(t)ri(t)Mt))y 2 <\m\- 
Therefore 

b b 
L(a) = J |Q«( t )(Q!,Q!)[ 1/2 dt < J \r{t)\dt = \r{b)\=L(a). 


The equality holds if and only if r(t) is monotone and (g a (t)Ti(t), T](t)) = 0. In this case 
the velocity of a satisfies the equation a(t) = f(t)d(exp p ) w u\u(t) = f(t)g a u\£(a(t)). From 
the other hand the extremal 7(5) = exp p (r(t)u(to)) , s = r(t), satisfies the equation 7(5) = 
d(exp p ) r M u r to \u(to) = s<7 7 ( s )£(7(s))- Since a(t) and 7(5) satisfy the same equation and have 
the same initial point, we conclude that a is a reparameterization of the extremal 7. □ 

We have noticed that a general piecewise smooth horizontal curve does not have a unique 
cotangent lift. If the curve is an extremal then there is a special cotangent lift, the one 
that satisfies the Hamilton-Jacobi equation. In the case of the two step bracket generating 
distributions it is possible to find a canonical cotangent lift. The condition for this is formulated 
in the following lemma. 

Lemma 4.15. Assume the strong bracket generating hypothesis. Let x(t) be any Lipshitz 
horizontal curve. Then there exists such a cotangent lift (x(t),£(i)) that a cotangent vector 

1 dg pq 

u>j = £j H — -£, P £,q is orthogonal to T J (^,v(x)) for any v £ S% at a.e. t so that pairing 

2 ux 3 

(ujj, r j (£, v(x))) = 0. This cotangent lift does not depend on the coordinate system and is 
called a canonical cotangent lift. 

Proof. Let (x(t),r](t)) be any cotangent lift and let . . . , y( n_m ) be a basis of sections of 
null-space S 1 - over a neighborhood of the curve. Since £ belongs in general to T*/S^r x S% 
then we can write 

n—m 

Z(t)= V (t)+v(t)=n(t) + Y J *k(t)v^(t). 

k=l 

Then 

& + 2 Q x j ^ = + V 3 + 2~dxT^ P + Vp) ^ q + Vq) 

1 dgP q ( 1 dgP q 

Recall that for w € S^r 

duo 1 ' duj ' 

Here we used that g pq £p-K-^- = ^q~^~ = Wj- Notice that T fc (£, w) = T k (r], w). We had already 
shown that T 3 (£, w) transforms as a tangent vector. Show now that 
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+ 



1 dgwjx) 

2 dxi 



Then 



£, p t; q transforms as a cotangent vector. Consider the transformation laws 



dx k 



£j and g- 



,pq 



QyP Qyl 

dx r dx s ' 



nPQ . 



^ + ^^rUq = 9 p %^ + 7 : 



2 <V' 



0& , lSflP 9 
p dy q 



2 % 



rs dy p dy q dx l d 



dx k \ dx m ^ 1 d 
dx r dx s dy p S1 dx m \ dyi J dy q 2 dx l 



dx l dx k dx n 



<9x r <9x s / Sy- 3 <9y p ch/ 9 



rs <9y p <9y 9 <9x' <9 2 x fc <9x m rs dy p dy q dx l dx k d^k dx m 



dx r dx s dyP dx m dyi dy q 
^ 1 dg rs dy p dy q dx l dx k dx n ^ 1 rs 
2 dx l dx r dx s dyi dy p dy q 2 



dx r dx s dyP dyi dx m dy q 
d 2 y p dy q dy p d 2 y q \ dx l dx k dx 



+ 



dx l dx r dx s dx r dx l dx s J dyi dy p dy q 



Qx k dt^jc 1 3g rs dx k 

The second and the third terms here are equal to g rs ^r^-^ — 6- and - - , Tr^£r£s respectively, 

oyi ox s 2 ox K oyi 

which gives in the whole the transformation of a covariant vector: 



dx k 
dyi 



a& idsr 

1 ^ r dx* + 2dx k ^ s 



dx 

w 



The rest of the terms give in sum 



2„k 



d 2 x 



i 



, , s d 2 y p dx 1 dx k 
dx s dyi r 2 dx l dx s dyi dy p 



2 dx k 



1 d 2 y p dx 1 dx k 



8 2 



x 



dx s dy- 



. d 2 y p dx 1 dx k 
dx l dx s dyi dy p 



&t.k = g rs ir£,k 



2 J dx l dx s dyi dy p ^ rsk 

d 2 x k d 2 y p dx 1 dx k 
dx s dyi dx l dx s dyi dy p 







since -tt^S* = 0, where 5 k is a Kronecker symbol. 

Now we see that the orthogonality condition is of the form 



ldg 



pq 



1 ij + 2^ x J r]pr]q ) + \^ + o 



ldg 



pq 



2 dxi 



■VpVq 



■T k (rj,w) = 0. 



As r(£, •) is injective the converse matrix (T k (r],w)) 1 exists. Therefore the linear system of 
n — m equations in n — m variables afc(i) 



k=l 

is uniquely solvable. 



1 finPI 



k=i 



VpVq) 



□ 



5. Differential of the exponential map 

As it was mentioned, the exponential mapping exp p is not a diffeomorphism at the origin, 
but as in the case of sub-Riemannian geometry there is a hope that it is a local diffeomorphism 
at some points. The main result can be stated that the exponential map exp p (u) is a local 
diffeomorphism if u is neither a null vector no an annihilator. We consider only the case 
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of 2-step bracket generating distribution. First, let us set out the Taylor expansion for k-th 
component of exp p (u), where p is fixed at the origin of the coordinates and u E T*M: 

N 1 

(5.1) exp p (n) fc = ^ ^- pr u Pl ...u Pr + 0(\u\ N+1 ), 

r=l 

where 7,^ "' Pr is symmetric in indexes p\,...p r and will be computed later, \u\ is any Euclidean 
norm on T*M. Notice that ex.p p (tu) = x(t), where — solution of the system (14,3ft 

with x(0) = 0, £(0) = u. Then at the origin 



We count for some value of t 

r+l 



Qk P1 ... Pr 

where £ Pi (t) denotes the absence of £ Pi (t)- Now, using (|4.3p and changing indexes, we get 
j\r+l Q kpi-Pr 

jt) xk{t) = {- J ^« t ))-9 qpr+r {m 



dtj V 

, f) a PrPr + l 



1 fjri^ryr + L \ 



2 'W dxi 

Therefore, 
(5.2) 

Q kpi...p r 

l%r PT+1 (x) = sym(pi, . . . tPr+1 ) • [9^(x)^-^(x) - ffi-M( x )W£l( x ) 

here sym{p\, . . . ,p r +i) means that we symmetrize the indexes pi, ■ ■ ■ ,p r +i- Setting r = in 
the previous formula we get 

(5-3) 7?) = = 9 kp , 

dx k 

since 7^ = x k (0) and — - equals to 1 if and only if k = q and zero otherwise. 



Analogously, observe that for r = 1 in (|5.2 



= lymfan) • (V 2 ^ - \9 kqd ~£f) = ~r kPlP2 - 

It is rather hard to calculate a general term, but it will be sufficient for us to look into the 
view of 7(3) . 

Now (|5.ip receives the following form 

k N 1 

exp p (u) fc = 7(f) 1 %i + Yl V\\r)"' Pru P± ■ ■ ■ u Pr + °(\ U \ N ) 



Hi) 1 k,dff 



P\P1 



r=2 
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and differentiating it, we obtain 

N 1 

(5.5) de Wp (u) k = t$ + £ _- 7 ^-^ 1Va ...u Pr + 0(\u\ N ) 



(r-1)! 



r=2 

N 



(r-1)! 

r=2 v ; 



More precisely, 



1 

2 

1 kjpq n 



dexp P (u) k * = g kj (0) + 7 k i) 2 u P2 + -^ 2P3 u P2 u P3 + 0(\u\ 3 ) 



g kj {0) _ F k jPup + ^ UpUq + (|„ 



Since we assumed 2-step bracket generating hypothesis, choose coordinates near p so that p 
is an origin and 



gi k (0) = / £ > 



EjP k 



where V is a m x m unit matrix and fi^/- 7 is a m x m matrix with v negative unities on 
the diagonal and m — v positive unities, which can be also written as follows: g^ k (0) = £j8 3 k , 
where 5 k is a Kronecker symbol and 

f-1, if 1 

£j = < 1, if v < j ^ m, 
I 0, if m < j ^ n. 

Denote with a, b the indexes responsible for elements standing in rows or columns with 
numbers 1, . . . , m, and a, (3 — for m + 1, . . . , n respectively. Then dex.p p (u) is a n x n matrix 
of the following form 

W kj = (A ab B a P 
~ \C ab D af3 

with 

A ab = £ jab + 0(\u\), 

(5.6) B a P = -T a ^u p + 0(\u\ 2 ), 

D^ = \^ q u p u q + 0(\uf). 

Since 7^ p = due to the special choice of g k i there are no terms of order 2 in D a/3 . The 
following proposition is an easy computation on determinant. 

Lemma 5.1. detW^(u) = detP^(u) + 0(|u| 2 ( n ~ m ) +1 ), where WW(u) is obtained from W by 
discarding the error terms containing 0(\u\ l ), i = 1, 2, 3 and detTy(u) is homogeneous of 
degree 2(n — m) in u. 
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To estimate the determinant of W(u) we need some more calculations. From (I5.2p and (15.4 

we get 



7 ( f 9 (*) 



a l{2) 



sym(/3, p,q) • ( g jg (x)-^-(x) - ^ ) 3 {x)^—{x) 



sym((3,p,q) ■ -g jq {x 



8xi (*) + r^)^(.) 



3 V 1 J 



aO + r apj (x) ^— (x) + r Q « ^- ( x i 



<9x-? 



9 n {x)^—(x) - g^(x)^—{x) - g^{x) 



dxi 



dxi 
dxo 



(x) 



Setting here x = 0, we get that the first and the last terms in the last sum are zero, because 
gi k (0) = for j, k > m. Hence, for p, q ^ m 



(5.7) 



aBab 
T(3) 

aBpq 



3 \ <9x-? <9x.? <9x b a <9x a 



and for p, q > m 7^ P9 = since g jk (0) = for j, k > m. Let us calculate the involved terms 
in (|577l) . 



x=0 



8x b 



x=0 



x=0 



,ak 



Q g aj 

8x k 

1 

2 



8g a i 
8x k 



,jk d 9 c 



8x k 



,ak 



.dg 00 



3k 



8g° 



x=0 

8g a(i 



8x a 



dxi 



,,ak 



8x K 8x K ax K 

\ ( Qg<*k dg Ba ak 8 2 g l3a 8gP k 8g' 



x=0 

cm 



1 8g af3 
2 £a ~8x^' 



2 V 8x b 8x k ^ 8x b 8x k 8x b 8x k 

Q2 g aa Q g ak Q g a8 Q2 g a8 , 



+ 9 



8k _ 



8x b x k 8x b 8x k 



9 



8x b x k 



x=0 



owing to g ak = g@ k = and 
Take a null-section v(x), then 



lfdg ak 8gP a 8g? k 8g aa 8 2 g af3 
2 \ 8x b 8x k 8x b 8x k a 8x a 8x b 

8g al3 



8x k 



by Lemma 13.21 Now we simplify the form of 



8 2 g a ' 3 
8x a 8x b 



dvi 



-Vk- 



8xP 8xP 

Thus, making use of Lemma 13.21 and differentiating both parts, we obtain 

9 ( g ik dv k\ = dg^dv^ + njk 8 2 v k 
8x q \ 8x p J 8x q 8x p 



8 ( 8g 



8x1 \ Q X P 



Vk 



Q2gjk 

dxi 8xP 



Vk + 



8x p 8x q ' 

Qgjk Q Vk 

8xP 'dx q ' 



From here 



8 2 gi k (x] 
8x q 8x p 



Vk{X) 



8g3 k {x) 8v k (x) 8g jk (x) 8v k (x) jk 8 2 v k (x) 

9 J (x)- 



8x p 8x q 



8x q 8x p 



8x p 8x q 



22 ANNA KOROLKO, IRINA MARKINA 

Taking inner product with another null-section, we get 

d 2 g jk (x) dg ]k {x) dv k {x) dg 3k (x) dv k (x) 

— V k (x)Wj(xj = Wj(x) Wj(X) 

dxidxP v ; n ; dxP dxi n ' dxi dxP n ' 

since o- 3 (x)— — - — Wj(x) = by virtue of Lemma I3T21 

y v 1 dxPdxi J 1 — ik 

Set x = and, since </ fc (0)|^(0) = ej|^(0) = -ei^-(0)v k , then 



dx r dx r dx 

d 2 gi k _dg^dvk gg^ dv k 

dx a dx bVkW '' dx a dx bW ^ dx b dx a ^ 



_ dg jk dg km _ dgi k dg km 

£k dx a dx b mW ^ £k dx b dx a Vm 



3- 



Take v k (0) = 5^, Wj(0) = 5 J a and get 



Q2 g af3 Q g . 



dx a dx b dx a dx b dx b dx a 
Therefore, 

dT af3a _ 1 ( dg aj dgP a dg? j dg aa dg aj dg^_ dg aj dgM 

dx b 2 \ dx b dxi Qx b dxi a 3 Qx a dx b a 3 dx b dx a 

Substituting calculated terms in (15, 7j) 

a/ 3ab _ If / dg a i dg aa \dgP b ( dg a i dg ab \dg^ a 

7 ( 3 ) ~ 6[{ £a ~dx^ +£j ~dxT) ~dxT~ \ Eb ~d^ + £] ~dxT) ~dxT 

dg aj Q g pa g g /3j dg aa g g aj Q g f3j Q g aj Q g f3j 



' dx b dxi Qx b dxi EaE] Q x a dx b £a£j dx b dx a 
dg aj dgP b dgM dg ab dg aj dg^i dg aj dg^ j 

dx a dxi dx a dxi £b£ i dx b dx a £b£ ^ dx a dx b 

To simplify this let us introduce the following notations 

E = £b ^ Ub ^ F " = ^ Ub - 



Then 



ctBab 
7(3) u a u b 



i UejE^E^ - •2/-."'/-/' - -jl-Jl-)' + E^If) 

6 (M* - EPj ) {( F ? - £ i Eaj ) - 3 M* + 



6 

- £j B jP B ja + 2e j B jP C aj . 
3 



Thus, we have the form of the matrix W i 



, £iP b B a $ 
W kj = ~ 1 ~ ~ ~ ~ 
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therefore, 



from which we obtain 



; - ° i - h pt i f. 

-C ab EjPP) \ ^B^B^ 



detVFl = \det-B j/3 B ja \ 
i i 3 



From here we have the homogeneity of det W(u) of degree 2(n — m) in u, since the matrix B 
is represented by the mapping F(u, •) : S 1 - — > S, where S 1 - is (n — m)-dimensional. 

Lemma 5.2. Let us assume 2-step bracket generating hypothesis for the ss-manifold M , and 
let u € T*M. Then for every u with (gu,u) ^ there exists 5 > such that 

|detM(u)| > 8\(gu,u)\( n -' m K 

Proof. By Theorem 13.11 the mapping T(u, •) : S 1 - — ► S is injective for every nonzero u with 
(gu,u) ^ 0. From the other hand, -B- 7 " is a matrix for — T(u, •) by (|5.6p and, hence, the 
matrix for B^B^ a is the matrix for r(u, •)* r • T(u, •), which is invertible by injectivity of 
r(n, •). Therefore, detP^(u) ^ if (gu,u) ^ and the statement of the lemma holds due to a 
homogeneity argument. □ 

Remark 5.1. Lemma 15.21 can be reformulated in the following way: detiy(u) ^ if and only 
if gu is a 2-step bracket generator. 

Theorem 5.2. i) If gu is a 2-step bracket generator, then there exists 5 > such that exp p (tu) 
is a local diffeomorphism for any < t < 5. 

ii) Assuming 2-step bracket generating hypothesis, there exists 5 > depending continuously 
on p such that exp p (u) is a local diffeomorphism for u near p = and (gu, u) ^ 0. 

Proof. The assertion i) follows from Lemma 15.21 and the Remark 15.11 By Lemma |5. II 

det W(u) ^ detW(u) - C\u\ 2{n - m)+1 

( \(gu,u)\\ {n ~ m) 

for small u. Thus, ii) holds for \u\ ^ 8 I — j— — j • □ 



6. Quaternion ss-manifold 

In the present chapter we find the parametric equations of extremals for a group furnished 
with the sub-semi-Riemannian metric of the index 2 described earlier in Example 2 at Section 
2. The Hamiltonian function H (£, 9, x, z) has the following form 

1 1 

H = - il + £ + £f) + -{x 2 x A 9 l e 2 + (x 2 x 3 + xix 4 )9i0 3 - x^O^) 

+\{0l-0l-9l)(-x\-xl + xl + xl) 
(6.1) +2^1 (-^26 + + x 46 - Z3&0 + ^0 2 {x^i + x 3 & + X 2 & + ^l&O 

+^3(^36 - £46 + a?i6 - ^2^4)- 
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The corresponding Hamiltonian system is 

-6 - 2^201 + \x462 + 5X3(93, 
-6 + \xi0i + \x 3 e 2 - \x40 3 , 





oh 


X\ = 


Wi 




dH 


X2 = 


W2 




oh 


X 3 = 






dH 


X4 = 





(6-2) 

£3 + \x±0i + \x 2 2 + ^16*3, 



U - \x 3 0i + - hx 2 03, 

a ?4 2 2 2 

dH 

h = -QQ- = \{{x 2 x A - xix 3 )9 2 + (x 2 x 3 + xix 4 )0 3 ) + \0\{-x\ - x\ + x\ + x|) 

+5 (-X26 + + x 4 £ 3 - x 3 6), 
dH 

i2 = "SZ - = \( x 2Xa ~ xix 3 )9i - \0 2 {-x\ -x\ + x\ + x\) + §(x 4 £i + X362 + X 2 ?3 + xi6), 
at/2 

^ = = 5(^2X3 + xix A )6i - \9 3 {-x\ -x\ + x\ + xf) + i(x 3 £i - x 4 6 + xi£ 3 - x 2 6), 

£l = - = -5(-X3^1^2 + X4M3) + - 6 2 2 - 61) - ±60! - ^ 4 2 - ^3^3, 

6 = = 4(x 4 0i0 2 + X3M3) + \ x 2{el - e\ - e\) + \^e x - \&e 2 + \&o z , 

ox 2 
d H 

£3 = ~k = -\(x 2 e x e 3 - x x e x e 2 ) - \x 3 (o\ - el - el) + ±£ 4 0i - \&e 2 - ^3, 
0x3 

dH 

£4 = = -i(x 2 0102 + XM) - 1*4(0? " 2 2 " 3 2 ) " 5^01 " 2^2 + ^203, 



dx4 
dH 

dz\ 
dH 

dz 2 

dH 

dz 3 



0, 
0, 
0. 



We observe that 01,02,03 are constants. Let us remind that the projection of a solution of 
the Hamiltonian system onto (x, z)-space is called extremal. In order to find it we will reduce 
the Hamiltonian system to the system containing only (xi, x 2 , x 3 , X4, zi, z 2 , z 3 ) coordinates. If 
we express £i> • • • >?4 from the first 4 equations and substitute them in the equations of the 
Hamiltonian system, then we obtain 

£1 = ^(x20i - x 4 2 - x 3 3 ), 
£2 = ^(-xi0i - X 3 6»2 + x 4 3 ), 

£3 = ^(X401 + X 2 02 + Xl0 3 ), 
£4 = ^(-X 3 0l + Xi0 2 - X 2 03)- 
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Differentiating first 4 equations and substituting £1, . . . ,£4 there, we get 

Xi = -x 2 1 + x 4 6> 2 + x 3 6> 3 , 

X 2 = X\Q\ + X3O2 - X4^3; 
X 3 = X 4 6>i + X 2 #2 + £l#3j 
X 4 = -±3^1 + £l#2 - ^2^3 

or 

/ -e x 



(6.3) 



P3 
02 





12 = 9i 

We are looking for the solution x\ = x\(t), . .. ,x& = x 4 (t), t E [—00, +00], satisfying xi(0) 
0, . . . , 24(0) = and £i(0) =£§,..., i 4 (0) = x 4 . The eigenvalues of the matrix 



X2 
X3 
\ X 4 J 



( 



A 



\ 











\ 



/ 



are Ai = a, A2 = —a, A3 = a, and A4 = —a, where a = \k\ + a 
k = 62 — 1O3. The associated eigenvectors are 

v\ = (ia\k\,a\k\,ak,iak), 

V2 = (— ia\k\, a\k\, —ak, iak), 

V3 = (ia\k\, —a\k\, —ak,iak), 

v 4 = (— ia\k\, —a\k\,ak, iak), 



\k\ — i9\ and k = 02 + ids, 



where \k\ = + 0§. Notice that the matrix A is skew-symmetric with respect to our 
nondegenerate metric Q with index 2) in the sense that Q ■ {Ax){y) = —Qx ■ Ay. This extends 
the idea of sub-Riemannian case, which was considered in [2], where the matrix Q was just a 
unit matrix and A was skew-symmetric in the usual sense. Also it carries on the sub-Lorentzian 
case, where A was skew-symmetric with respect to sub-Lorentzian metric Q [11]. 
The solution of the system (|6.3p is of the form 



where 



(6.4) 



xi(t) 
x 2 (t) 

X3{t) 
Xi(t) 

c\ 

(■■2 
<:■■ 

C4 



c 2 ae~ at + c 3 ae at 



c 4 ae ~ at ), 



i\k\( Cl ae at 
\k\(aae at + c 2 ae~ at - c 3 ae m - c 4 ae~ m ), 
c\ake at - c 2 ake~ at - c 3 a£;e a * + c 4 ake~ at , 
i(c\ake at + c 2 ake~ at + C3ake m + c 4 ake~ M ) 



1 



4iak\k\ 
1 

Aiak\k\ 
1 

4iak\k\ 
1 

Uak\k\ 



{k(xi + ij%) + \k\(x% + ia%)), 



-fc(i$ 



IX 



^ + \k\(4- i 4)), 



(k(x»-ix») + \k\(xi-i4)), 
{-k(x1 + ix 2 ) + \k\(xl + ix 3 )). 
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Therefore, the x-coordinates of the extremals have a form 

xi{t) = i\k\( Cl e at + c 2 e~ at + c 3 e m + c A e~ m ) - i\k\(a + c 2 + c 3 + c 4 ), 

(6.5) x 2 (t) = \k\{ Cl e at - c 2 e~ at - c 3 e m + c^~ m ) - \k\{a - c 2 - c 3 + c 4 ), 
xz{t) = cike at + c 2 ke~ at - c 3 ke m - ake~ m - (c\k + c 2 k - c 3 k - 04k), 
x A (t) = i{c\ke at - c 2 ke~ at + c 3 ke at - c^ke~ at ) - i(cik - c 2 k + c 3 k - c 4 k). 

From the horizontality conditions 

1, 

Zl = -{+X 2 Xl - XlX 2 + X4X3 - X3X4), 

I, . 

Z 2 = -(-X4X1 - X 3 X 2 + X 2 X 3 + X1X4), 

1, 

Z 3 = -(-X3X1 + X 4 X 2 + X1X3 - X 2 X 4 ) 

we can find the vertical components 

(6.6) z x (t) = 2i\k\ 2 (-2( Cl c 2 a - c 3 c A a)t + dc 2 {e at - e~ at ) - c 3 c 4 (e M - e~ M )), 
z 2 (t) = 20 2 |A;|(-2(c 1 c 2 a + C3C4a)t + ciC2(e a *-e- a *)+C3C4(e a -e- a )) 

+ 28 1 6 3 (c 1 c 3 e mt + c 2 c 4 e~ 2|fc|t - C1C3 - c 2 c 4 ) 

+ 2ie 3 \k\(c lC3 e at + c 2 c 4 e- at - Cl c 3 e m - c 2 c 4 e- m ), 

z 3 {t) = 26 3 \k\(-2( Cl c 2 a + c 3 c 4 a)t + Cl c 2 (e at - e~ at ) + c 3 c 4 (e m - e~ m )) 

- 26>i6» 2 (cic 3 e 2|fc| * + c 2 c 4 e~ 2|fc|t - C1C3 - c 2 c 4 ) 

- 2id 2 \k\{ Cl c 3 e at + c 2 c 4 e- at - Cl c 3 e m - c 2 c 4 e- M ). 

The constants of integration a, i = 1, ... ,4 are given by (16. 4D through the initial velocity. We 
would like to calculate the homogeneous norm of an element (x(t), z(t)) given by 

||(x,z)|| 4 = (-x 2 -x 2 +x 2 + x^) 2 + z 2 + z 2 + z 2 . 

We have 

(6.7) ||x(t)|| 2 = (_x 2 -x 2 + x 2 + x 2 )(t) 

= 8\k\\2 Cl c 2 + 2c 3 C4 - c lC2 (e at + e~ at ) - c 3 c A {e at + e""*)) 

= — 32|/c| 2 (cic 2 sinh 2 + c 3 c 4 sinh 2 y ). 

Let us introduce the notation 101 = ^(x? + ^c 2 )> w 2 = £4 + ix-j in order to simplify the 
calculations. Then 

Cl ° 2 = ~ l6a 2 \k\ 2 C^ 2 ' 2 " ^ + 2ilm ( w ^) 



and 



We see that c\c 2 = c 3 c 4 and (|6.7p takes the form 



C3C4 = - 16 _ 2|fcj2 (\w 2 \ 2 - I 2 - 2ilm(w 1 w 2 )) . 



\x(t)\\ 2 = -64|/c| 2 Re(cic 2 sinh 2 y), a= |fc|+»0i. 



We also need the values 
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C2C4 = -wmw lW2 - wi12 ' 



and 



ClC 2 C 3 C 4 = - l62|Q|4|fc|4 

Then 

{z\ + zl + zl){t) = 4|A;| 4 (l6t 2 ciC2C 3 c 4 aa - 8cic 2 c 3 c 4 at{e at - e' at ) - 8cic 2 c 3 c 4 at(e s - e~ m ) 
+4cic 2 C3C4(e a * - e~ at )(e m - e~ m ) + c?c|(2e 2|fc| * - e 2at - e 2m ) + clcl^ 2 ^ 1 - e~ 2at - e~ 2m ) 

-ic lC2 c 3 c 4 + 2c lC2 c 3 C4(e 2iei * + e" 2 ^)) 

+48 2 \k\ 2 (cic 2C3 C4(4 - 2(e 2 l fc l' + e^')) - 2( Cl c 2 e 2 ^ + c 2 c\er mt ) 



+ ( c 2 c 2 e 4|^ + c 2 c 2 e -4|fc| t) + (c 2 c 2 + c 



2 C 4/ 



+8^i|fc| 3 (c 2 C 2 ( e ( a+2 l fe l)* - e ("+ 2 l fe D* - e a * + e m ) + clcKe'^ 2 ^ - e^ 2 ^ - e~ at + e - st ) 

+C l C 2 C Z C A {(e {a - 2 ^ t + e -( a -2|fc|)*) _ ( e (S-2|fc|)t + e -(a-2\k\)tj _ ^at + g -at) + ^ot + e "^))^ . 

To simplify the last expression we notice that 

16|fe| 4 ciC2C3C 4 (l6t 2 aa - 8t(ae at - ae- at + ae m - ae" at ) + (e 2 l fc l* + e" 2 ^ - e 2i(>lt - e" 2 * 1 *)) 
= 64cic 2 c 3 C4|A;| 4 (t 2 (|A;| 2 + # 2 ) - 2i(|fc| sinh(|fc|t) cos(0i*) + 0\ cosh(|lfe|t) sin(^t)) 

+ sinh 2 (|A;|t) +sin 2 (^t)) 

= 64cic 2 c 3 c 4 |A;| 4 ((|fc|* - sinh(|fc|t) cos(<M)) 2 + {9 x t - cosh(|fe|i) sin(<M)) 2 

and 

4|fe| 4 (c 2 c 2 (2e 2 l fc l* - e 2at - e 2m ) + c 2 C 2 (2 e - 2 l fe l' - e~ 2at - e~ 2m ) 

-4cic 2 c 3 c 4 + 2cic 2 c 3 c 4 (e 2iei * + e^ 1 ')) 
= 4|fc| 4 (4sin 2 (0 1 t)(c 2 c|e 2 l fc l i + c^e" 2 ^ 4 ) - 8ci c 2 c 3 c 4 sin 2 (^)) 
= 16[A;| 4 sin 2 (e 1 t)(c 1 c 3 el fc l t - c 2 c 4 e-^ t ) 2 , 

4^ 2 |A:| 2 ( ClC2C 3C 4 (4 - 2(e 2 l fc l* + e" 2 ^)) - 2(c l c 2 e 2 ^ t + c^e" 2 !*"*) 
+(c 2 c 2 e 4 l^ + c 2 c 2 e-^) + (c 2 c 2 + c 2 c 2 )) 
= 4# 2 |£f( - 8 Cl c 2C3 c 4 sinh 2 + c 2 c 2 (l - e 2 ^ 1 ) 2 + c 2 c 2 (l - e " 2 l fc l*) 2 ) 

= 16ef|ife| 2 siiili 2 (|A;|t)(ciC8e |fc| *-C2C4e- |fc| *) 2 , 



and 
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and 

8^ 1 |fc| 3 ( C 2 C 2( e ( a + 2 l fc l)* - e (W\W _ e at + e ») + c 2 c 2 (e -(a+2|fc|) t _ e -(a+2\k\)t _ g -ai + g-at) 
+CiC 2 C3C4((e( a - 2 l fc l)* + e -(«-2|fc|)i) _ ( e (S-2|fc[)i + e -(3-2|fc|)ij _ ( g at + g -at) + ^at + £ -at^ 

= 8^i|A;| 3 (4i sinh(|fc|t) sin(0it)(c? c|e 2|/c| * + c^e -21 * 1 *) - 8icic 2 c 3 c 4 sinh(|fe|t) sin(0it)) 
= -320i|fc| 3 smh(|fc|t)sm(0it)(cic 3 el fc l* -c 2 c 4 e~l fc l') 2 . 
Then joining last expressions we get 

(z\ + zl + zl)(t) = 16|Af(cic 3 e |fc| * - c 2 c 4 e- |fc| ') 2 (#i sinh(|fc|t) - \k\ sin^t)) 2 
+64|A;| 4 ci c 2 c 3 c 4 ((\k\t - smh(\k\t) cos((M)) 2 + {6 it - coeh(|fc|t) sin(6M)) 2 ) . 
The letter calculations we can summarize in the following 

Theorem 6.1. Let Q be ss-manifold described in Example 2 of Section 2. The normal ex- 
tremals starting from the identity of the group Q are given by the solution of the Hamiltonian 
system (|6.2p by the parametric equations ()6.5[) and (|6.6p . The constant of integrations ct, 
k = 1,2,3,4, are related to the initial velocity (x®, x 2 , x 3 , x 4 ) by ()6.4j) . where k = 62 + i03 and 
a = \k\ + %B\ are the first integrals of the Hamiltonian system (|6.2|) . 
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